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Chains of rotational tori and filamentary structures close to high 
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This paper discusses phase space structures encountered in the neighborhood of periodic orbits 
with high order multiplicity in a 3D autonomous Hamiltonian system with a potential of galactic 
type. We consider 4D spaces of section and we use the method of color and rotation [Patsis and 
Zachilas 1994] in order to visualize them. As examples we use the case of two orbits, one 2- 
periodic and one 7-periodic. We investigate the structure of multiple tori around them in the 
4D surface of section and in addition we study the orbital behavior in the neighborhood of 
the corresponding simple unstable periodic orbits. By considering initially a few consequents in 
the neighborhood of the orbits in both cases we find a structure in the space of section, which 
is in direct correspondence with what is observed in a resonance zone of a 2D autonomous 
Hamiltonian system. However, in our 3D case we have instead of stability islands rotational tori, 
while the chaotic zone connecting the points of the unstable periodic orbit is replaced by filaments 
extending in 4D following a smooth color variation. For more intersections, the consequents of 
the orbit which started in the neighborhood of the unstable periodic orbit, diffuse in phase space 
and form a cloud that occupies a large volume surrounding the region containing the rotational 
tori. In this cloud the colors of the points are mixed. The same structures have been observed 
in the neighborhood of all m-periodic orbits we have examined in the system. This indicates a 
generic behavior. 
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1. Introduction 

Recently, Katsanikas and Patsis [2011] (hereafter 
KPll) studied the structure of the phase space in 
the neighborhood of simple periodic orbits in a 3D 
autonomous Hamiltonian system of galactic type. In 
the present paper we extend this work and we inves- 
tigate the orbital behavior in the neighborhood of 
periodic orbits with high order multiplicity. Espe- 
cially, we study the orbital behavior in the neighbor- 
hood of a stable 2-periodic and a stable 7-periodic 
orbit and in the neighborhood of the accompanying 
simple unstable periodic orbits. 

In Cartesian coordinates {x, x,y,y, z, z), if we 
consider our surface of section to be defined by 
y = 0, a "2-periodic" orbit is one that closes af- 
ter 2 intersections with the surface of section when 
y > 0. This is a periodic orbit of multiplicity 2. In 
the same way, in general, an "m-periodic" orbit is a 
periodic orbit of multiplicity m. Here we study two 
cases with m = 2 and 7. For the visualization of the 
4D surfaces of section we use the method of color 
and rotation [Patsis and Zachilas 1994]. With this 
method we plot the consequents in a 3D subspace 
of the 4D surface of section and every consequent is 
colored according to its location in the 4th dimen- 
sion. For a description of this method, the meaning 
of viewing angles etc., see KPll. 

Our Hamiltonian is of the form 

H{x,y,z,x,y,z) = 
^{x"^ + f + z^) + ^{x,y,z) 

-l^li^' + y') (1) 

where ^{x, y, z) is the potential we used in our 
applications, i.e.: 

^{x,y,z) = 

GMi 

"(x2 + ^ + [ai + (f^ + 6f)V2]2)i/2 " 

GM2 

(x2 + |^ + [a2 + (4 + 6i)V2]2)l/2 
ya yj, 

(2) 

This potential is the same as in KPll, 
i.e. a triaxial double Miyamoto disk rotating 
around its short axis z with angular velocity 
Qb=60 kms~^kpc~^. 

In our units, distance R=l corresponds 
to 1 kpc. The velocity unit corresponds to 



209.64 km/sec. For the Jacobi constant Ej=l cor- 
responds to 43950 {km/secf'. For the rest of the 
parameters we have used the following values: ai = 
kpc, hi = 0.495 kpc, Mi = 2.05 x 10^° M©, 02 = 
7.258 kpc, 62 = 0.520 kpc, M2 = 25.47 x 
lO"*^*^ Mq, Qa = 1.2, qb = 0.9. The parameters qa,qb 
determine the geometry of the disks, while a, b are 
scaling factors. 

The trajectories are calculated numerically by 
using a 4*^^ order Runge-Kutta scheme, with a con- 
stant time step that secures a 13 digits precision. A 
typical integration for 10^ consequents with a Core 

2 Duo CPU/2. 2GHz computer takes about 3.6 min- 
utes of real time. 

The calculation of the linear stability of a peri- 
odic orbit is based on the method of Broucke [1969] 
and Hadjidemetriou [1975]. By this method we cal- 
culate the stability indices bl, b2 and the quantity 
A (following the notation of Contopoulos and Mag- 
nenat [1985]). Depending on the values of the sta- 
bility indices and that of A, a periodic orbit can 
be stable (S), simple unstable (U), double unstable 
(DU) and complex unstable (A). For definitions see 
Contopoulos and Magnenat [1985]. For a general- 
ization of this method to systems with higher than 

3 degrees of freedom see Skokos [2001]. 

2. Spaces of Section 

2.1. Orbits 

The stable and simple unstable orbits of the 2- 
periodic case belong to the 3D families of 2-periodic 
orbits "s" (an initially stable family) and "u" (an 
initially simple unstable family). They are bifurca- 
tions of a planar 3/1 family on the equatorial plane. 
The morphology of the orbits of the families s and 
u is depicted in Figs. [Hand [2] respectively. 

We call the 3D 7-periodic families si (an ini- 
tially stable family) and ul (an initially simple un- 
stable family). They are bifurcations of the family 
xlvl [Skokos et al 2002a,b] that is associated with 
the vertical resonance 2/1 in our galactic system. In 
Figs. [3l m we show the morphology of the orbits of 
the families si and ul respectively. 

2.2. spaces of section in the 
neighborhood of the 2-periodic 
orbits 

2.2.1. Phase space structure close to s 

We consider first the stable 2-periodic orbit s at 
Ej=-4.33035. We 
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Fig. 2. Projections of a simple unstable orbit of the 2-periodic family u at Ej= —4.33035. 



trace it at the initial conditions {xq, zq^xq^ zq) = 
(1.0867145, 0.2555081, -0.56256937, 0.049693274). 
In the neighborhood of s we observe two tori, when 
we perturb its initial condition xq, by Ax belong- 
ing to the intervals 10~^ < Ax < 4 x 10~^ and 
— 10"'^ < Ax < — 5 X 10~^. All these perturbations 
of the initial conditions are isoenergetic displace- 
ments on the surfaces of section. For example if we 
take Ax = 4 x 10~^ we observe two tori in Fig. [5] 
around the two points (black dots) of the stable 2- 
periodic orbit. On the two tori we observe a smooth 
color succession. This means that the succession of 
the colors of the consequents on a given structure 
follows the succession of the colors on the color bar 
at the right side of the figure. The consequents are 
colored according to the value of the coordinate, 
that is not used in the 3D spatial projection. In all 
diagrams of this paper the color is given according 
to the z coordinate. The color values in the color 
bars are normalized to the [0,1] interval correspond- 
ing to [zmin, Zmax]- Foi the details of the method 
the reader may refer to KPll. The smoothness of 
the color variation can be checked by looking at the 
color bar, on the right of Fig. [5l where we see that 
the color variation is between neighboring shades. 



The colors on the upper left torus vary from red to 
orange to yellow and to green, while the colors on 
the other torus vary from green to light blue to blue. 
This means that the fourth value of the consequents 
has a smooth distribution in the 4th dimension. If 
we perturb the initial conditions with larger values 
of perturbations, for example Ax = —6 x 10~^, we 
observe again two tori with smooth color variation 
for 1000 consequents as in the previous case. How- 
ever, if we continue the integration of the orbits we 
observe that the consequents start to deviate from 
two tori and they form a cloud in the 3D projection 
of the 4D surface of section around them (Fig. [6]). 
In the cloud the color is mixed and this means that 
the points are far away in the 4h dimension. The dy- 
namical behavior according to which an orbit stays 
close to an invariant torus for some time and then 
diffuses in phase space is typical of sticky orbits 
[Contopoulos and Harsoula 2008] . In Table 1 we give 
the values and the direction of the perturbations for 
which we observe tori. Just beyond this interval we 
encounter sticky orbits as the one in Fig. [6l 
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Fig. 3. Projections of a stable orbit of the 7-periodic family si at Ej= —4.622377. 
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Fig. 4. Projections of a simple unstable orbit of the 7-periodic family ul at Ej= —4.622377. 



2.2.2. Phase space structure close to u 

Close to s, at Ej= —4.33035, we trace also the 
simple unstable 2-periodic orbit u (Fig. [5]). Our 
code finds it at initial conditions {x(),zq,xq,zq) = 
0.67649546, 0.16719438, 0.78985004, 0.41351876). In 
its neighborhood we have found two types of dy- 
namical behavior. We can see the first type if we 
perturb for example the initial conditions by Ax = 
10"^. In this case we observe in Fig. [7] a "ribbon" 
that connects the points of the simple unstable 2- 
periodic orbit. If we apply the method of color and 
rotation we observe in Fig. [8] a smooth color vari- 
ation from red to blue. This means that the "rib- 
bon" is a 4D structure. After 5300 intersections we 
observe in Fig. O that the consequents leave the 
"ribbon" and they scatter in the phase space. This 
orbit is also sticky in the sense that remains close 
to a given phase space structure for some time and 
then diffuses in the 4D space of section. 

The second type of dynamical behavior can be 
seen in the neighborhood of the simple unstable 2- 
periodic orbit for Ej= —4.33035. If we add for exam- 
ple a perturbation in the z-direction Az = 2 x 10^'^ 
in the initial conditions of u we observe a filamen- 
tary surface in the 3D projection of the 4D surface 
of section (Fig. llOp . This filamentary structure con- 
nects the points of u and forms four loops. The two 
of them surround the two tori that are around the 



two points of s. In the regions, close to the points 
of u we have two self-intersections of the filamen- 
tary structure. In Fig. [10] we observe two more self- 
intersections indicated with arrows and the forma- 
tion of 4 loops. The loops surround the periodic 
points s and two more points representing a periodic 
orbit symmetric to s with respect to the equatorial 
plane z = 0, while the arrows represent an orbit 
symmetric to u. The dynamics in the neighborhood 
of this second set of periodic orbits is similar to the 
one close to s and u. Now we apply the method of 
color and rotation in order to study the distribution 
of the consequents in the 4th dimension (Fig. [TT]) . 
We observe, that we have also in this case a smooth 
color variation from red to blue and this means that 
the consequents are on the filamentary structure in 
the 4D space of section. For more than 4700 inter- 
sections the points deviate from this structure and 
form a cloud of points around it. Therefore in this 
case we encounter again the phenomenon of stick- 
iness. Table 2 gives the range of the perturbation 
of the initial conditions for which we find the "rib- 
bons" or the "filaments" connecting the points of 
u, for Ej = —4.33035. For larger values of the per- 
turbations we find clouds of points as the one we 
presented in Fig. [9l 
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2.3. 4D spaces of section in the 
neighborhood of 7-periodic 
orbits 

The structures observed close to the 2-periodic or- 
bits, have been found in the neighborhood of ev- 
ery m-periodic orbit we have studied in our system. 
The results are qualitatively similar, but the fila- 
mentary structure that joins the points of the sim- 
ple unstable orbits becomes complex as m increases. 
Below we give one more example of the phase space 
structure close to a high multiplicity orbit, this time 
around a 7-periodic one. 



2.3.1. Phase space structure close to si 

We apply a perturbation both in the x-direction and 
z-direction in the initial conditions of the stable 7- 
periodic orbit of the family si. Ax = 1.6 x 10~^ 
and Az = -5 X 10"^ for Ej= -4.622377. At this 
Jacobi constant we find si at initial conditions 
{xo,zo,xo,zo) = (0.25992717,0.29991487,0,0). The 
perturbed orbit has been integrated for 10^ con- 
sequents. In the (x, X, z) 3D projection of the 4D 
surface of section we observe seven tori surrounding 
the points with the initial conditions of si (Fig. I12p . 
In this projection we observe, that two of the tori 
intersect each other. To check their position in the 
4th dimension we apply the method of color and 
rotation (Fig. [T2]) . By taking into account the val- 
ues of all tori in the 4th dimension, z, and scaling 
the colors according to [imim •2maa:]=[— 0.617, 0.617], 
we find that to each small torus correspond one 
or mostly two primary colors. We observe that in 
the case of the two tori that intersect each other 
in the 3D projection (x, x, z) we have the meeting 
of different colors at their intersection. This means 
that this intersection does not exist in the 4D space 
but only in the 3D projection, as expected. Next 
we apply the method of color and rotation only to 
one torus, which we call Tl and which is indicated 
with an arrow in Fig. [T2j In Fig. [13] we observe a 
smooth color variation on its surface, which obvi- 
ously corresponds to a different z range than that 
in Fig. [T2| and we see a succession of colors from red 
to orange, to yellow, to green, to light blue, to blue 
and finally to violet. Comparing this result with this 
found in KPll we realize that we have a morphol- 
ogy typical for rotational tori [Vrahatis et al 1997, 
KPll]. The only difference is that we do not ob- 
serve in the present case the transition of the color 
sequence from the external to the internal surface of 




Fig. 5. The 4D surface of section for Ej=-4.33035 (3000 
consequents) in the neighborhood of the stable 2-periodic or- 
bit s. The initial conditions of s (black dots) are perturbed by 
Aa; = 4 X 10~^. We use the (x, x, z) projection to depict the 
consequents and the i value (—0.253 < i < 0.708) to color 
them. We give also the initial conditions of the associated 
simple unstable 2-periodic orbit u with red dots. Our point 
of view in spherical coordinates is {9,<j)) — (22.5°, 0°). 



the torus, as observed in cases around stable simple 
periodic orbits in KPll (compare with figure 11 in 
KPll). 

If we perturb the initial conditions of the peri- 
odic orbit of si by Ax = 10~^ we see in Fig.ll4lseven 
tori without any intersection of these tori in the 3D 
subspace. These tori have also a smooth color vari- 
ation (Fig. [T4l) . In Table 3 we give the range of the 
perturbation, for which we find tori in the neigh- 
borhood of si. 



2.3.2. Phase space structure close to ul 

Now we investigate the properties of orbits 
in the phase space in the neighborhood of 
the simple unstable periodic orbit ul. The or- 
bit we study is at Ej = —4.622377 and is 
found with initial conditions {xq, zq, xq, zq) = 
(0.27756296, , 0.29935776, -0.084579242, 0.16982513). 
Firstly we add a perturbation in the initial condi- 
tions in the x-direction Ax = 10~^. In Fig. [15] we 
see that the consequents that are depicted with red 
color double bow in the 3D projection of the 4D sur- 
face of section that surrounds the seven tori that are 
around the points of si. If we color the consequents 
according to the value of their fourth dimension we 
see in Fig. [16] that we have a smooth color variation 
from red to violet but at the two intersections A and 
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Fig. 6. 1500 consequents of the sticky orbit we found in the 
4D surface of section for Ej=— 4.33035. Color is given to the 
consequents according to their z value (—1.99 < i < 1.982). 
The orbit stays initially close to the region of the two tori 
(indicated with arrows) and then diffuses in phase space. We 
have chosen a {6, <j)) = (22.5°, 22.5°) point of view for a better 
inspection of the regions where the consequents stay initially 
concentrated and the way they diffuse. 




-1.5 



Fig. 7. The 3D subspace {x,x,z) of the 4D surface of sec- 
tion for Ej=— 4.33035 in the neighborhood of a and u. We 
apply a perturbation Ax = 10 to their initial conditions 
and we consider 5000 consequents. The initial conditions of a 
are indicated with "+" in green color. We give also the initial 
conditions of u with x symbols in black. Our point of view 
in spherical coordinates is {9,<j)) — (26°, 298°). Around s a 
double torus structure is formed, while the perturbed orbit 
in the neighborhood of u forms a filamentary structure. 

B we have the meeting of different colors (e.g. blue 
with orange) . This means that we have different val- 
ues of the fourth dimension at the two intersections 
and these intersections are only projection effects 
and not true intersections in the 4D space of sec- 
tion. The consequents depart from the structure we 
give in Fig. \T6\ after 25000 intersections. 

If we perturb the same orbit by Az = —3.6 x 
10-"^ at Ej= -4.622377. We observe in Fig.ITl that 
the consequents (with red) form a filamentary struc- 




Fig. 8. The orbital behavior close to u for a deviation 
Ax = 10^"* from its initial conditions, at Ej = —4.33035 
(5000 intersections) in the 4D surface of section. We use 
the {x, X, z) space for plotting the points and the i value 
(—0.637 < i < 0.638) to color them. Our point of view in 
spherical coordinates is given by (6, (f)) = (180°, 22.5°). 




Fig. 9. A (x, X, z) 3D projection of the 4D surface of sec- 
tion close to the periodic orbits s and u (x symbol) for 
Ej=— 4.33035. The "ribbon" we observe in red is the one 
given in 4D in Fig. [S] while the blue structures are the tori 
around s. This time we integrate the orbit of Fig. [S]to get 
7500 (instead of 5000) consequents. We observe a cloud of 
points, that surrounds both the tori and the "ribbon". Our 
point of view in spherical coordinates is (9, <j)) — (26°, 240°). 

ture that connects the points of the simple unstable 
7-periodic orbit and surrounds the seven tori of the 
stable 7-periodic orbit. In Fig. [17] with the numbers 
1 until 7 we indicate the 7 points of ul in the {x, x, z) 
subspace. At these regions the filaments, that have 
been formed by evolving the orbit in time, cross 
each other in the 3D projection of the 4D surface 
of section (Fig. [T7|l . We observe, that besides the 
crossings at the numbered regions, where we have 
the initial conditions of ul, we have 7 more crossings 
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Fig. 10. The same projection as in Fig.O but with an orbit 
close to u perturbed this time in the ^-direction by Aa: = 
2 X 10""^ (4500 consequents). This time the "ribbon" has 
four self-intersections in the 3D subspace. We indicate s with 
green "+" symbols, u with "x" in red, while the arrows point 
to the two additional self-intersections of the "ribbon". Our 
point of view is {9, 0) = (26°, 298°). 
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Fig. 11. The 4D representation of the "ribbon" of Fig. [TOl 
for 4000 intersections. The consequents are colored according 
to their i value (—0.651 < i < 0.651). The branches of the 
"ribbon" that intersect themselves, have the same color at the 
intersections, meaning that they are true intersetions in the 
4D space. Our point of view is given by (6, 0) = (22.5°, 22.5°). 



of the filamentary structure by itself and the forma- 
tion of 7 new loops, like those surrounding the tori 
around the points of si. We have again here, like in 
the 2-periodic orbit in section 2.2.2, the presence of 
the symmetric family, with respect to the equatorial 
plane. 

In Fig. [18] we observe a smooth color varia- 
tion along the consequents that build the filaments. 
Starting from 1 we have a color succession from 
green to light blue, then to blue at 2, after that 
to light blue, then to green at 3, etc. At the regions, 




0.26 



Fig. 12. The 4D surface of section for Ej=-4.622377 (6500 
consequents) in the neighborhood of si. The initial conditions 
of si are given with black dots. The orbit of the tori is found 
by perturbing the si initial conditions (see text) by Aa; = 
1.6 X 10~^ and Az = — 5 x 10~^. The {x,x,z) projection is 
used to depict the consequents and the i value (—0.617 < 
z < 0.617) to color them. We give also the initial conditions 
of the associated simple unstable 7-periodic orbit ul with red 
dots. Our point of view in spherical coordinates is (6, (j)) = 
(0°, 180°). We indicate with an arrow the Tl torus, which is 
discussed in the text. 



0.245 



0.306 



0.28. 




0.06 



Fig. 13. Torus Tl in the 4D surface of section. Color is given 
to the consequents according to their z value (—0.092 < i < 
0.092). We observe a smooth color variation on its surface. 
Our point of view in spherical coordinates is given by (9, (f>) — 
(22.5°, 45°). 



where we have the self crossing of the filaments, i.e. 
close to the points 1, 2,. ..7, we observe that the re- 
gions are characterized by just one color only. For 
example at 1 the shade is green, at 2 blue etc. This 
means that the fourth coordinate z of the conse- 
quents at these regions has this time the same value 
and the intersections in the 3D projections are real 
intersections in the 4D space. We also observe that 
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Fig. 14. Tori in the 4D surface of section for Ej=-4.622377 
(10^ consequents) in the neighborhood of the stable 7- 
periodic orbit si for Aa; — 10^^ . Color is given to the conse- 
quents according to their i value (—0.592 < i < 0.592). Our 
point of view in spherical coordinates is given by {9, (j>) ~ 
(22.5°, 45°). 



Fig. 16. The 4D surfaces of section in the neighborhood of 
ul for the same orbit close to it depicted in Fig. 1151 We con- 
sider 6000 intersections. We use the {x,x,z) space for plot- 
ting the points and the z value (—0.708 < i < 0.693) to color 
them. Different colors at the intersections (at "A" and "B") 
indicate that they are not true in 4D. Our point of view is 
(e,<^) = (22.5°,45°). 




Fig. 15. The 3D {x, x, z) projection of the 4D surface of 
section for Ej =— 4.622377 at the neighborhood of si and 
ul. Both si, ul initial conditions have been perturbed by 
Ax — 10~^ to get a set of seven tori and a "ribbon" struc- 
ture respectively. Green "-|-" correspond to si and black 
"x" to ul. Our point of view in spherical coordinates is 
(e,,^) = (28°,48°). 

the color evolves along both branches that depart 
from, or arrive at the points of the unstable periodic 
orbit. Both branches show the same color evolution 
between two successive crossings. 

The smooth succession of the colors along the 
filaments is observed as long as the consequents par- 
ticipate in the filamentary structure. However, if we 
consider more than about 6500 consequents in the 



3 




Fig. 17. The 3D {x,x,z) projection of the 4D surface of 
section for Ej=— 4.622377 at the neighborhoods of si and 
ul for a perturbation of their initial conditions by Az = 
—3.6 X 10^^. The initial conditions of si are indicated with 
"x" and green color and those of ul with "x" black symbols. 
These points are numbered from 1 to 7. Our point of view in 
spherical coordinates is (0,(f>) = (26°, 192°). 

neighborhood of the ul orbit, the points diffuse, oc- 
cupying a large volume in the phase space and fi- 
nally they form a cloud in the 3D projections that 
surrounds the structure we observe in Fig. [181 The 
dimensions of this cloud becomes about 15 times 
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Fig. 18. The 4D representation of the orbit close to ul, 
which we give in Fig. 1171 for 6000 intersections. We use 
the {x,x,z) space for plotting the points and the z value 
(-0.726 < z < 0.709) to color them. The points 1 to 7 in- 
dicate the ul initial conditions. The common colors of the 
branches of the filamentary structures that meet at the in- 
tersections show, that they are true intersections in 4D. Our 
point of view in spherical coordinates is given by {6, (f)) = 
(180", 22.5°). 

larger in the x direction than the structure depicted 
in Figs. [T7] and [121 This is determined by the space 
in which the particles are allowed to move for this 
Ej. The cloud of points can be seen in Fig. [191 The 
red configuration in the central region of the di- 
agram is the structure we observe in Fig. [181 By 
giving colors to the consequents according to their 
distribution in the 4th dimension we have realized 
that in the cloud the distribution of points in the 
4th dimension is mixed (Fig. [20l) . In Table 4 we 
summarize the range of perturbation for which the 
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Fig. 19. The 3D {x,z,x) projection showing a cloud of 
points that surrounds the orbits close to si and ul presented 
in Figs. [17] and 1181 This time we integrate the same orbit 
in the neighborhood of ul as in the previous figures, but we 
consider 10000 instead of 6000 consequents. The red struc- 
ture in the center of the cloud includes everything plotted in 
Figs. [T71 Our point of view in spherical coordinates is given 
by (e,0) = (37°,17°). 




Fig. 20. The cloud of Fig. [19] in the 4D space of section. We 
plot 7500 consequents. We use the {x,x,z) space for plotting 
the points and the z value (—2.29 < i < 2.33) to color them. 
We observe mixing of colors. Our point of view is (6, <j)) = 
(180°, 22.5°). 

consequents remain on the filamentary structure. 

3. Lyapunov Characteristic 
Numbers 

In this paper we study two kinds of orbits. The first 
kind is represented by tori in the 4D surface of sec- 
tion in the neighborhood of stable 2-periodic and 
7-periodic orbits. The second kind is represented 
initially by structures confined in phase space (dou- 
ble bow, ribbon, "filamentary" structures etc.) and 
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Table 1. The intervals of the perturbations of the initial conditions of the s orbit at Ej — —4.33035, for which 
we find 4D tori around s. 



Aa; 


Ax 


Az 


Ai 


from 10"* to 4 X 10"^ 


from 10"* to 3 X 10"^ 


from 10"* to 10"'' 


from 10"* to 6 X 10"^ 


from -10"* to -5 x 10"'^ 


from -10"* to -3 x 10""" 


from -10"* to -10"'' 


from -10"* to -1.1 X 10"^ 



Table 2. The range of perturbations in the initial conditions of u, that give ribbon or filamentary structure in the 4D 
surfaces of section for Ej — —4.33035. The first two rows correspond to the ribbons and the two last to the filaments. 



Aa; 


Ai 


Az 


Ai 


from 10"* to 3 X 10"-^ 


from 10"* to 2 X 10"^ 


from 10"* to lO""* 


from 10"* to 3 X lO""* 


from -10"* to -10"-^ 


from -10"* to -3 x 10"^ 


from -10"* to -4 x 10"-* 


from -10"* to -9 x 10"* 


from 4 X lO""* to 6 x 10"^ 


from 3 X 10"^ to 6 x 10"^ 


from 2 X lO""* to 10""^ 


from 4 X lO""* to 4 x 10"^ 


from -2 X lO""* to -5 x 10"^ 




from -5 X lO""* to -2 x 10"^ 


from -10"-^ to -3 x 10"^ 



then by clouds in the 4D surfaces of section. These 
orbits are located in the neighborhood of simple un- 
stable 2-periodic and 7-periodic orbits. In this sec- 
tion we calculate the "finite time" Lyapunov Char- 
acteristic Number (LCN) for these two types of or- 
bits. 

The "finite time" Lyapunov Characteristic 
Number is defined as: 



1 



LCN{t) = - In 



m 



where ■^(to) and (,{t) are the distances between two 
points of two nearby orbits at times t = and t 
respectively (see e.g. Skokos [2010]) 

Firstly, for the first kind of orbits, we com- 
puted the maximal Lyapunov Characteristic Num- 
ber (mLCN). For example, for the orbit that is rep- 
resented in the 4D surface of section in Fig. [T^ the 
value of LCN(t) decreases like 1/t and tends to zero 
as we can see in Fig. [2TJ 

Then, we calculated the "finite time" Lyapunov 
Characteristic Number {LCN{t)) for the second 
type of orbits, for example for the orbit in Figs. [18] 
and Fig. [T9l In Fig. [18] we have for the first 6500 
consequents a filamentary structure. During this pe- 
riod the LCN(t) of the orbit decreases to a value 
1.4 X 10"^ (Fig. [22li). Beyond that point the orbit 
is represented by a cloud of points in the 4D surface 
of section. During this phase, the LCN(t) fluctuates 
as time increases and flnally increases and tends to 
level off around 1.67 x 10"^ (Fig. [22b). 

4. Conclusions 

In this paper we studied the phase space structure 
in the neighborhood of stable and simple unstable 
periodic orbits in a rotating 3D galactic potential. 
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Fig. 21. The LCN(t) for the orbit given in the 4D surface 
of section in Fig. 1141 The axes are in logarithmic scale. 
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Fig. 22. a)The evolution of LCN(t) for the orbit that is 
presented in Figs. [18] and Fig. \T9\ (b) The part of Fig. [22b 
for 70000 <t< 115000. The axes are in logarithmic scale. 



We presented the dynamical behavior of two sets 
of stable-simple unstable orbits. The flrst was in 
the neighborhood of periodic orbits of multiplicity 
m = 2, and the second in the neighborhood of pe- 
riodic orbits of multiplicity m = 7. For less than 
6000 consequents we observed a direct correspon- 
dence between the standard conflguration, which we 
encounter in a resonance zone in 2D autonomous 
Hamiltonian systems, and the structure of phase 
space in our 4D spaces of section, namely a succes- 
sion of elliptic and hyperbolic points. This means 
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Table 3. The range of perturbations in the initial conditions of si, that give tori or tori with intersection in the 3D pro- 
jections of the 4D surfaces of section for Ej = —4.622377. The first two rows correspond to the tori and the two last to 
the tori with intersection. 



Aa; 


Ai 


Ax 


Ai 


from 10~* to 1.3 x 10"'' 


from 10"* to 3 X 10"^ 


from 10"* to 1.6 x 10"^ 


from 10"* to 10"^ 


from -10"" to -1.1 X 10"^ 


from -10"* to -3 x 10"^ 


from -10"* to -5 x 10""^ 


from -10"* to -10"^ 


from 1.4 X 10"^ to 1.8 x 10"^ 


from 4 X 10"^ to 6 x 10"^ 






from -1.2 X 10"'' to -1.7 x 10"^ 


from -4 X 10"^ to -6 x 10"^ 







Table 4. The range of perturbations in the initial conditions of wl, that give double bow or filamentary structure in the 4D 
surfaces of section for Ej = —4.622377. The first two rows correspond to the double bows and the two last to the filaments. 



Aa; 


Ai 


Az 


Ai 


from lO"'^ to 7 X 10"* 


from 10"* to 3 x lO""" 


from 10"* to 2 X 10"* 


from 10"* to 6 X lO""* 


from -10"* to -9 x 10"* 


from -10"* to -4 x lO"'^ 


from -10"* to -2 x 10"* 


from -10"* to -9 x lO"'^ 


from 8 X 10"* to 1.6 x 10"^ 


from 4 X lO""* to 2 x 10"^ 


from 3 X 10"* to 4 x lO""* 


from 7 X lO""* to 10"^ 


from -10^'' to -10"^ 


from -5 X 10"'' to -8 x 10"^ 


from -3 X 10~* to -5 x lO"-" 


from -10"^ to -10"^ 



that also in the case we study we observed a suc- 
cession of stable and simple unstable points in the 
4D spaces of section. For larger integration times 
the consequents diffuse and occupy a larger volume 
of the phase space. We found similar behavior for all 
cases of m-periodic orbits we studied in this system. 

The following are the main conclusions from our 
work: 

• In the neighborhood of the stable m-periodic or- 
bit we found m tori surrounding its initial condi- 
tions. We found smooth color variation along these 
tori in two ways. Along a particular torus and also 
along all m tori, considering them as one object. 
This depends on the scale of the fourth coordinate, 
which gives the colors to the consequents. These tori 
arc rotational tori in the terminology introduced by 
Vrahatis et al. [1997]. Similar structures can be seen 
in the paper by Martinet & Magnenat [1981]. How- 
ever, we did not observe in the present study the 
color transition from the external to the internal 
side of the torus on the individual tori as in some 
cases in KPll. 

• Integrating an orbit close to a simple unstable m- 
periodic orbit to obtain a few thousands of con- 
sequents, we found that they form a filamentary 
structure with smooth color variation in the 4th 
dimension. The filaments that are formed in this 
way connect the points of the simple unstable m- 
periodic orbit and surround the seven tori around 
the points of the stable m-periodic orbit. The fila- 
mentary structures appear either as "ribbons" or as 
"bows" in the 4D spaces of section. 

• In the regions close to the points of the simple un- 



stable periodic orbit two branches of a filamentary 
structure meet and their consequents have the same 
color. This shows, that we have at these points self- 
intersections of the structures in the 4D space of 
section. 

• The above described situation is a direct extrapola- 
tion of the typical case of a 2D autonomous Hamil- 
tonian system with the chain of stability islands and 
the chaotic zone, which connects the points of the 
corresponding unstable periodic orbit. We note that 
in this 3D case the unstable orbit is simple unsta- 
ble and we have a smooth color variation along the 
filaments when wc consider the first few thousands 
of the intersections in the space of section. 

• We encountered cases, where the consequents that 
remained on the filaments and reinforced this struc- 
ture for a few thousand intersections, diffused later 
in the 4D space. The diffusion in the 4D space, was 
characterized by mixing of colors. 
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